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ABSTRACT 

We present integrated polarization properties of nearby spiral galaxies at 4.8 GHz, and models for 
the integrated polarization of spiral galaxy disks as a function of inclination. Spiral galaxies in our 
sample have observed integrated fractional polarization in the range < 1% to 17.6%. At inclinations 
less than 50 degrees, the fractional polarization depends mostly on the ratio of random to regular 
magnetic field strength. At higher inclinations, Faraday depolarization associated with the regular 
magnetic field becomes more important. The observed degree of polarization is lower (< 4%) for more 
luminous galaxies, in particular those with L4.8 > 2 x 10^^ W Hz"-'^. The polarization angle of the 
integrated emission is aligned with the apparent minor axis of the disk for galaxies without a bar. 
In our axially symmetric models, the polarization angle of the integrated emission is independent of 
wavelength. Simulated distributions of fractional polarization for randomly oriented spiral galaxies at 
4.8 GHz and 1.4 GHz are presented. We conclude that polarization measurements, e.g. with the SKA, 
of unresolved spiral galaxies allow statistical studies of the magnetic field in disk galaxies using large 
samples in the local universe and at high redshift. As these galaxies behave as idealized background 
sources without internal Faraday rotation, they can be used to detect large-scale magnetic fields in 
the intergalactic medium. 
Subject headings: galaxies: 
ization 



magnetic fields — galaxies: spiral — radio continuum: galaxies — polar- 



1. INTRODUCTION 

Large-scale regular magnetic fields have been detected 
by Faraday rotation of radi ation of polarized background 
sources m the Milky Way (|Han et al.lll999 l: iBrow n et all 
2003, 200 '^, and the nearest galaxies fStcpano v et al.| 
2008t lArshakianet all 12007^. such as M3l'(H an et al. 
, and the Large Magellanic Cloud (G aensler et al.l 
. The limiting factor for this work is the num- 
ber density of bright polarized background sources. A 
few galaxies in the Local Group (M31, LMC) have an- 
gular sizes large enough to probe their magnetic field 
structure by rotation measures of polarized background 
sources with current radio telescopes. Magnetic fields 
in galaxies outside the Local Group have been studied 
mainly with imaging polarimetry at radio wavelengths 
(Beck 2005 for a review). These observations indicate a 
large-scale aziniuthal magnetic field in the plane of the 
disk of spiral galaxies, normally oriented in a spiral pat- 
tern with a pitch angle similar to the optical spiral arms. 
High degrees of polarization have been observed in inter- 
arm regions of some galaxies, such as NGC 6946 (|Beckl 
[2007b). If a magnetic field is observed in a barred galaxy, 
it is aligned with the bar inside the regi on of the bar, wit h 
a spiral pattern outside the bar region ()Beck et al.ll2002l ). 
However, some barred galaxies show no detectable polar- 
ized emission or polarized emission along the bar. 

Statistical information on magnetic fields in galaxies 
is required to study the cosmic evolution of these mag- 



netic fields by observing galaxies over a range of redshift. 
Polarimetric imaging of a large sample of spiral galaxies 
in the local universe is impractical with current instru- 
mentation because of the high sensitivity and angular 
resolution required. On the other hand, deep polariza- 
tion surveys of the sky are underway or being planned 
that will have the sensitivity to detect spiral galaxies out 
to large distances, but not resolve these galaxies. The 
first deep wide-area polarimetry survey with the sensi- 
tivity to detect spiral galaxies is GALFACTS^ at 1.4 
GHz, 5(7 point source sensitivity of 400/iJy, and reso- 
lution 3^ SKA path finder telescopes such as ASKAP 
(jJohnston et al.ir2007[ ) aim to do all-sky polarimetry sur- 
veys at 1.4 GHz with an instantaneous band width of 
300 MHz with thousands of spectral channels and a 5g 
point source sensitivity around 50 /^Jy ( Johnston et al.l 
l2007t) . ASKAP will be able to detect a few thousand 
spiral galaxies at sufficient signal to noise ratio to detect 
polarization at the few percent level, but the expected 
angular resolution of 8" to ~ 1' will not resolve these 
galaxies. This raises the question what can be learned 
about magnetic fields from polarimetry of large samples 
of unresolved spiral galaxies. 

Observations of polarization from spiral galaxies at 
redshifts up to 1 or higher will require the Square Kilo- 
metre Array^ (SKA), and even then it will be diffi- 

1 http://www.ucalgary.ca/ras/GALFACTS 
^ http://www.skatelescope.org 
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cult to make resolved polarization maps of these high- 
redshift galaxies. The origin and evolution of cosmic 
magnetism is one of fi ve key science drivers for the SKA 
(|Gaensler et alj l2004j) . For the evolution of magnetic 
fields in spiral galaxies as a function of redshift, polar- 
ization of the integrated emission is most readily observ- 
able, and a comparison with unresolved galaxies at low 
redshift is relatively straightforward. 

The number density of faint polarized radio sources in 
the sky determines the unknown confusion limit in polar- 
ization, and is an important parameter for rotation mea- 
sure studies. Below 1 mJy total flux density, star forming 
galaxies become an increasingly important fraction of all 
radio sources. Knowledge of the polarization properties 
of these galaxies is required to make realistic predictions 
of the number of polarized sources at ^Jy levels in po- 
larized flux density. Spiral galaxies may become an im- 
portant fraction of the polarize d 1.4 GHz radio source 
population at /itJy flux densities (jStil et al.ll2007l ). 

This paper presents observations and models that show 
that unresolved spiral galaxies have detectable polar- 
ized emission, sometimes more than 10% polarized at 
4.8 GHz, and that statistical information on magnetic 
fields and internal Faraday rotation can be derived from 
integrated polarization measurements of spiral galaxies. 

2. OBSERVATIONS 

We use a compilation of archival data for nearby galax- 
ies imaged in polarized radio emission. Observations at 
4.8 GHz are most widely available, so this paper focuses 
at this frequency to maximize the sample size for a sin- 
gle frequency. Additional observations at 8.4 GHz are 
available for five galaxies. These are listed alongside the 
4.8 GHz data. Three sub samples are distinguished in 
this paper. The first sub sample is a set of 14 nearby 
"field" galaxies listed in Table [1] The second sub sam- 
ple is a set of 9 spiral g alaxies in the Virgo cluster from 
Wezgowiec et al.M 2007t ). Wezgowiec et al. (in prep.), and 
Chvzv et a l. (2007) listed in Tabled These galaxies ex- 
ist in a cluster environment, with a higher probability 
of interaction with other cluster members and the in- 
tracluster gas. The third sub sample of 20 barred spi- 
ral galaxies listed in Table [3] allows us to explore the 
integrated polarization of barred galaxies. One nearby 
galaxy (NGC 4631, classified as SB) was included in Ta- 
bled 

The nearest galaxies in the sample and the galax- 
ies in the Virgo cluster were observed with the 100m- 
Effclsberg telescope, sometimes in combination with the 
VLA. These galaxies are not affected by missing short 
spacings in the interferometer. The barred galaxies were 
not observed with the Effelsberg telescope, and the im- 
ages could in principle lack large scale structure resolved 
by the shortest baselines of the interferomet er. A com- 
pariso n of 6 cm flux densities derived by iBeck et al.l 
(|200^ with published integrated flux densities showed 
that NGC 1313 and NGC 3992 in Table[3]may be affected 
by missing short spacings, because published flux densi- 
ties listed in the NASA Extragalactic Database (NED) 
are substantially higher than the 6-cm flux density in 
IBeck et al.l()2002D . If no 6 cm flux density was published, 
we compared the 21 cm flux density correcte d for a spec- 
tral in dex 0.7 - 1 with the flux density in iBeck et aTl 
(|2002f ). Of the two barred galaxies with suspected miss- 



ing flux in the interferometer image, NGC 1313 showed 
no significant integrated polarization, while NGC 3992 
has a high degree of polarization of 14.9 ± 4.8 %. Most 
of the emission of NGC 3992 detected in the interferom- 
eter i mage originates outside the bar region (jBeck et al.1 
I2OOI ). 

In order to simulate unresolved observations from our 
maps of nearby galaxies we integrated the Stokes I, Q, 
and U images of these galaxies separately in ellipses. 
Their axial ratio was defined by the outer isophotes 
in the Stokes / image and the ellipses include the en- 
tire galaxy up to the noise level of the Stokes I map. 
The integrated Q and U flux densities were used to 
calculate the integrated polarized flux density Sp — 

\/(Q'^ + U'^), that was correc t ed fo r polarization bias 
following ISimmons fc Steward (|1985l ) to obtain the in- 
trinsic polarized flux density Sp^. The correction for 
polarization bias did not have a measurable effect for 
most galaxies because of the high signal to noise ra- 
tio of the integrated flux densities. The integrated 
fractional polarization of each galaxy was calculated as 
Ho = Sp_o/S, where S is the total flux density, includ- 
ing thermal and non-thermal emission. The polarization 
angle 6'poi — \ arctan(J7/Q) of the integrated emission is 
expected to be aligned with the apparent minor axis of 
the optical galaxy. The sensitivity of the observations 
is included in the errors for Hq through the noise in the 
maps. The error estimates also include the uncertainty in 
the zero level, that dominates the uncertainty for galax- 
ies with a large angular size such as M 31. It is customary 
to list the position angle of the optical major axis, and 
also to rotate the observed E-vectors by 90° to obtain 
B- vectors in the plane of the sky. We therefore define 
— &poi + 90° for comparison with published position 
angles of the optical major axis, ^opt- The results are 
listed in Tables [J [2 and [3 

Bandwidth depolarization is negligible in these data. 
For a rotation measure as high as 300 rad m~^, the po- 
larized intensity is decreased by less than 1%. The in- 
tegrated polarized flux densities are not affected by sig- 
nificant residual instrumental polarization. Small values 
for the integrated fractional polarization are the result 
of integrating significant flux in Stokes Q and U, with 
nearly equal positive and negative contributions. The 
uncertainty in the polarization angle of the integrated 
emission is typically less than 5 degrees. 

Some galaxies (M 31, M 81, NGC 4565 in Table [J 
NGC 4192, NGC 4302, NGC 4535 in Table[2l NGC 1300, 
NGC 1493, NGC 3953, NGC 3992 in Tabled have in- 
tegrated fractional polarization Hq > 7%. Five of these 
galaxies (M 31, NGC 4565, NGC 4192, NGC 1493, and 
NGC 3992) are more than 10% polarized after integra- 
tion over solid angle. 

Five of the Virgo galaxies have additional observations 
at 8.4 GHz. The fractional polarization at both frequen- 
cies agrees to within a few percent for these galaxies. 
Small differences can be attributed to statistical errors, 
mainly at 8.4 GHz where the signal to noise ratio of the 
emission is lower. The differences in Ho at 4.8 GHz and at 
8.4 GHz per galaxy are much smaller than the differences 
between galaxies. This confirms that the variation in Ho 
between galaxies represents intrinsic differences between 
the galaxies. The polarization angles at 4.8 GHz and 
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TABLE 1 

Integrated polarization of nearby spiral galaxies 



Source 


Type 


i 


V 




Sv{Vi Mpc)'' 


Ho 




^opt 


Tel." 






(°) 


(GHz) 


(mJy) 


(mJy) 


(%) 


(°) 


(°) 




IC 342 


SARfrsVH 


25 


4.8 


856 


82 


2.2 ± 0.1 


50 ± 1 


39 




M 31 




76 


4.8 


1863 


9 


17.6 ± 3.5 


68 ± 4 


38 




M 33 


SA(s)cd 


56 


4.8 


1539 


11 


1.4 ± 0.4 


77 ± 6 


23 


E 


M 51 


SAbc 


20 


4.8 


380 


358 


2.8 ± 0.3 


-46 ± 2 


-10 


E+V 


M 81 


SA(s)ab 


59 


4.8 


385 


41 


7.3 ± 0.3 


-36 ± 1 


-28 


E 


M 83 


SAB(s)c 


24 


4.8 


809 


111 


2.5 ± 0.4 


9 ± 4 


45 


E+V 


NGC 253 


SAB{s)c 


78 


4.8 


2707 


82 


1.9 ± 0.2 


51 ± 2 


52 


E+V 


NGC 891 


SA(s)b? sp 


88 


4.8 


286 


264 


1.6 ± 0.1 


6 ± 1 


23 


E 


NGC 3628 


SAb pec sp 


89 


4.8 


247 


111 


4.0 ± 0.1 


105 ± 1 


104 


E 


NGC 4565 


SA{s)b? sp 


86 


4.8 


54 


85 


10.5 ± 0.3 


-43 ± 1 


-44 


E 


NGC 4736 


{R)SAB(rs)ab 


35 


4.8 


125 


27 


2.4 ± 0.4 


-86 ± 3 


-65 


E+V 


NGC 5775 


SAbc 


81 


4.8 


94 


670 


0.5 ± 0.2 


-31 ± 9 


-35 


V 


NGC 5907 


Sc 


87 


4.8 


72 


87 


2.2 ± 0.3 


-6 ± 3 


-25 


E 


NGC 6946 


SAB(rs)cd 


38 


4.8 


457 


224 


1.5 ± 0.5 


-82 ± 7 


60 


E+V 



^ Flux density scaled to a distance of 10 Mpc.^ Imaging telescope: E — EfFelsberg 100m; E+V — EfFelsberg 100m 
combined with the Ver y Large Array ( VLA). References: IC 342: IGravc & BeclfJ (11988D, M 31: IBerkhui isen et al.l 
(|200|), M 33: Tabatab aei et al.l II2007I) . M 51: IFletcher er al. (2004), M 81: Beck et al. (1981), M 83: [Beck (2007a[) 
from Beck, Ehle & Sukumar (unpubl.), NGC 253: IHeesen et al. (20(18), N GC 891, NGC 3628, an d NGC 4565: 
iDumke & Krause (1998), NGC 4736: [C hvzv & Buta (2008), NGC 5775: ITiillmann et al.l (120001 ). NGC 5907: 
iDumke et al.i (i200a) , NGC 6946: IBeck fc Hoernea 11996,1 : .Bech (i200Tbl) 



TABLE 2 

Integrated polarization of Virgo galaxies 



Source 


Type 


i 


V 


5, 


Ho 






Tel." 






n 


(GHz) 


(mJy) 


(%) 


(°) 


(°) 




NGC 4192 


SABb 


78 


8.4 


16 


14.5 ± 1.0 


-28 ± 1 


-27 


E 








4.8 


38 


11.9 ± 0.5 


-21 ± 1 




E 


NGC 4254 


SA(s)c 


42 


4.8 


167 


1.4 ± 0.1 


-46 ± 1 


68 


E+V 


NGC 4302 


SAc 


90 


8.4 


22 


5.8 ± 1.0 


-13 ± 4 


-2 


E 








4.8 


36 


8.4 ± 0.7 


-8 ± 1 




E 


NGC 4303 


SABbc 


18 


8.4 


107 


1.5 ± 0.3 


-13 ± 4 


-28 


E 








4.8 


189 


0.7 ± 0.1 


-65 ± 4 




E 


NGC 4321 


SAB 


30 


8.4 


66 


2.3 ± 0.4 


-58 ± 4 


-60 


E 








4.8 


109 


1.0 ± 0.2 


-52 ± 4 




E 


NGC 4388 


SAb 


82 


4.8 


74 


1.1 ± 0.1 


-60 ± 3 


-89 


E 


NGC 4501 


SAb 


61 


4.8 


115 


1.0 ± 0.1 


-27 ± 3 


-42 


E 


NGC 4535 


SABc 


44 


8.4 


20 


6.7 ± 1.1 


17 + 3 





E 








4.8 


38 


8.0 ± 0.4 


15 + 1 




E 


NGC 4654 


SABc 


56 


4.8 


51 


3.3 ± 0.2 


-47 + 1 


-59 


E 



^ Imaging telescope: E — Effelsberg 100m; E+V — EfFelsberg 100m combined with the 
Very Large Array(VLA). References: NGC 4192, NGC4 302, NGC4303: unp ublished im- 
ages courtesy M. Wezgowiec and M. Ur banik, NGC 4254: IChvzv et all (I2007D . NGC 4388, 
NGC 4501, NGC 4535, and NGC 4654: IWezgowiec et al.l (I200W) . 



8.4 GHz are also mutually consistent. In low-inclination 
galaxies, the polarization angle of the integrated emis- 
sion is more affected by asymmetries because the effect 
of an axially symmetric component on the polarization 
vanishes mostly. This can explain the difference in po- 
larization angle for NGC 4303 {i = 18°). The position 
angle of the major axis can also be affected by ellipticity 
of the disk. 

The fractional polarization of any particular galaxy 
depends on the direction from which the galaxy is ob- 
served, as the projection of the magnetic field on the 
plane of the sky and the path length of the line of sight 
through the disk change. Any particular galaxy can only 
be observed from one direction, but the variation of in- 
tegrated polarization properties with inclination can be 
revealed for a sufficiently large sample of galaxies. The 
wide range in inclination (20° < i < 90°) in the sample 
allows an effective analysis of the integrated polariza- 



tion properties of spiral galaxies as a function of incli- 
nation. It should be noted that the sample is not ran- 
domly distributed in inclination. The galaxies in Ta- 
bles [U [21 and [3] were observed to study their magnetic 
field configuration through imaging polarimetry. The 
sample is biased towards low inclination (face-on) and 
high-inclination (edge-on) galaxies because these were 
preferred for imaging of the magnetic field. In partic- 
ular, the sub sample of barred spirals in Table [3] contains 
only galaxies with inclination i < 60° , with the exception 
of the added nearby galaxy NGC 4631. Any sample of 
spiral galaxies with inclination derived from the appar- 
ent axial ratio appears deficient in galaxies with i < 20° 
because the disks of most spiral galaxies are not exactly 
circular (iBinnev fc de Vaucouleurs Ill98l]: lGrosb0]||1985l : 
iLambas et al.lll992f) . This effect is also apparent in the 
present sample, but it does not affect our analysis. 
Figure [T] shows IIq function of inclination. The 
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TABLE 3 

Integrated polarization of barred galaxies 



Source 


TvDe 




1/ 


S,, 


5,^(10 Mpc) 


112 


0-a 
^ a 




n 

"opt 


Tel.'' 






\ ) 


(GHz) 


(mJy) 


(mJy) 


1%) 


(°'\ 




\ ) 




NGC 1097 


SBbc(rs) 


45 


4.8 


155 


397 


(0.7 ± 0.8) 






-45 


V 


NGC 1300 


SBb(s) 


35 


4.8 


10 


38 


7.6 ± 2.0 


89 ± 


6 


86 


V 


NGC 1313 


SBc(s) 


38 


4.8 


22 


4 


(0.7 ± 1.3) 






-10 


A 


NGC 1365 


SBb{s) 


40 


4.8 


205 


740 


(0.3 ± 0.2) 






40 


V 


NGC 1433 


SBb{s) 


27 


4.8 


2 


3 


(4.7 ± 3.1) 






17 


A 


NGC 1493 


SBc(rs) 


30 


4.8 


2 


3 


13.6 ± 3.5 


29 ± 


5 


C 


A 


NGC 1559 


SBc{s) 


55 


4.8 


105 


235 


1.9 ± 0.1 


-50 ± 


1 


-115 


A 


NGC 1672 


SBb(rs) 


39 


4.8 


105 


236 


1.5 ± 0.2 


-55 ± 


3 


-10 


A 


NGC 2336 


SBbc{r) 


59 


4.8 


4 


35 


(3.6 ± 3.5) 






-2 


V 


NGC 2442 


SBbc{rs) 


24 


4.8 


71 


182 


1.1 ± 0.5 


51 ± 


9 


40 


A 


NGC 3059 


SBc(s) 


27 


4.8 


29 


56 


(0.8 ± 0.5) 






C 


A 


NGC 3359 


SBc{s) 


55 


4.8 


13 


29 


(0.7 ± 1.4) 






-10 


V 


NGC 3953 


SBbc{r) 


61 


4.8 


6 


14 


8.5 ± 2.7 


50 ± 


6 


13 


V 


NGC 3992 


SBb(rs) 


59 


4.8 


4 


8 


14.9 ± 4.8 


25 ± 


7 


67 


V 


NGC 4535 


SBc{s) 


26 


4.8 


12 


30 


6.4 ± 1.4 


8 ± 


4 


28 


V 


NGC 4631 


SB(s)d 


85 


4.8 


476 


268 


1.1 ± 0.1 


16 ± 


2 


86 


E 


NGC 5068 


SBc{s) 


29 


4.8 


14 


7 


(1.8 ± 1.6) 






-20 


V 


NGC 5643 


SBc{s) 


30 


4.8 


59 


116 


1.0 ± 0.2 


-68 ± 


4 


C 


A 


NGC 7479 


SBbc(s) 


45 


4.8 


32 


365 


3.9 ± 0.6 


-9 ± 


3 


25 


V 


NGC 7552 


SBbc{s) 


31 


4.8 


110 


485 


1.2 ± 0.2 


-87 ± 


4 


1 


A 



^ Galaxies with Ho less than 2 times the formal error listed in parentheses; these galaxies are not included 
in figures'^ Imaging telescope: A — Autralia Telescope Compact Array (ATCA), E — Effelsberg 100 m, V 
^ VLA.*^ Position angle of major axis undetermined for this low-inclination galaxy. References: NGC 4631: 
Krause et al. (unpubl.) but see |Krause| j200^ and lGolla & Hummell jl994f) ; all other galaxies from lBeck et all 




30 60 90 30 60 90 

Inclination (degrees) Inclination (degrees) 

Fig. 1. — (a) Integrated fractional polarization at 4.8 GHz of nearby spiral galaxies from Table [T] (filled circles) and Virgo cluster spirals 
from Table [2] (triangles) , with model curves described in Section |3] The model curves from top to bottom correspond with models 1-9 in 
Table [5l in that order, (b) The same as (a) for the detected barred galaxies from Table |3] (squares). The model curves in (b) are shown for 
reference only. 
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Fig. 2. — (a) Correlation between 9opt and 0b for galaxies in Table [T] (filled circles) and the Virgo galaxies in Table [2] (triangles) . The 
symbol size increases linearly with sin(i). (b) The same as (a) for the barred galaxies in Table [51 (squares) . 
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Fig. 3. — Integrated fractional polarization as a function of distance (a) and monochromatic 4.8 GHz luminosity (b) for the entire 
sample. The symbols are the same as in Figure|2] The symbols for the Virgo galaxies NGC 4321 and NGC 4501 overlap completely in (a), 
and partially in (b). 
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curves in Figure [T] represent models that will be discussed 
in Section [3l The formal errors on Hq are generally small 
because of the high signal to noise ratio in the integrated 
flux densities. Some galaxies with large angular size have 
larger errors because of uncertainties in the zero level of 
the maps. The large spread in XIq, 1% £ Hq < 15%, 
seen for each subsample reflects intrinsic differences. The 
highest values of Hq in Figure [1^ are found in the incli- 
nation range from 50° to 80°. Some low- inclination 
galaxies with high Hq are found among the barred galax- 
ies in Figure [Jd. The model curves in Figure [Jd are 
shown for reference only, as they do not include a bar. 

Figure [2k shows the correlation between the position 
angle of the major axis, 0opt, and for galaxies in Ta- 
ble [1] and Table [2l If polarization of the integrated emis- 
sion is the result of projection of the azimuthal magnetic 
field in the disk on the plane of the sky, a correlation is 
expected between 9b and the line of nodes of the disk 
intersecting the plane of the sky. The latter is deter- 
mined by the position angle of the major axis, 9b- The 
most deviating galaxies in Figure [5^ are NGC 6946 and 
NGC 4254, with inchnation i = 38° and i = 42° respec- 
tively. Large-scale departures from axial symmetry are 
most likely to be noticed for galaxies with a small incli- 
nation. The position angles used in Figure [2] are also less 
well defined. At lower inclinations, the position angle of 
the major axis is more difficult to measure because the 
isophotes are nearly circular. Besides measurement un- 
certainties, the position angle of the major axis is not a 
good estimator for the position angle of the line of nodes 
of the disk if the inclination approaches 26°, the incli- 
nation for which a circular disk would have an apparent 
axial ratio of 0.9, similar to the mean intrinsic e lliptic - 
ity of spiral galaxy disks found bv lLambas et al.l (|1992f ). 
Consequently, the uncertainty associated 0opt is substan- 
tial but hard to quantify for galaxies with inclination less 
than ~ 40°. For galaxies with higher inclinations, 6'opt 
is estimated to be accurate to within ~ 5°. The uncer- 
tainty in 9b is typically a few degrees, but some addi- 
tional scatter is expected because of differences in struc- 
ture between galaxies. The total scatter in Figure [2^ for 
galaxies with inclination > 60° suggests that the scatter 
introduced by structural differences is limited to ±15°. 

Figure [2)3 shows the same relation for the detected 
barred galaxies in Table [3] The scatter between 
and appears substantially larger than in Figure [2p,. 
This may be a result of magnetic fields aligned with the 
bar, where the orientation of the bar adds a parame- 
ter to the problem. In this paper we focus mainly on 
galaxies without (strong) bars, while noting that unre- 
solved barred galaxies can in fact also be polarized radio 
sources. 

Figure [3] shows the integrated fractional polarization 
Ho as a function of distance, and monochromatic lumi- 
nosity ^4.8- There should be no correlation between in- 
tegrated polarization and distance, and none is found in 
the present sample. Figure [Sb shows a correlation be- 
tween Hp and L4.8. There are no galaxies in the sample 
with L4.8 > 2 X 10^1 W IIz~i that have Hq > 4%. Most 
galaxies with L4.8 < 5 x 10^° W Hz"^ do have Hq > 4%. 
As shown in Figure [TJ a galaxy can have low IIq because 
of strong Faraday depolarization, or because of its incli- 
nation. There is no correlation between inclination and 
luminosity in the sample. The lower integrated polariza- 



tion of luminous galaxies in Figure [SJd suggests stronger 
Faraday depolarization in these galaxies, or a lower de- 
gree of uniformity of the magnetic field as will be dis- 
cussed in Section |47T] 

3. MODELS OF THE INTEGRATED POLARIZATION 
3.1. Polarized intensity in the disk 

The principal cause of net polarization of unresolved 
spiral galaxies is the regular azimuthal magnetic field in 
the disk projected on the sky along the apparent major 
axis. Faraday rotation and depolarization are affected by 
the path length through the disk, the line-of-sight com- 
ponent of the regular magnetic field, and fluctuations 
in thermal electron density and magnetic field strength. 
A first analysis of the integrated polarization of spiral 
galaxies is best served with the simplest possible model 
that includes these factors. We consider axially sym- 
metric models with a regular azimuthal magnetic field 
in the disk, and an isotropic random component of the 
magnetic field. 

Using this simple geometry, the integrated polariza- 
tion of a spiral galaxy is modeled by considering each 
solid angle element of the disk as a slab with a speci- 
fied line-of-sight component of the regular magnetic field 
that depends on inclination and location in the disk. 
ISokoloff et al.l (|l998) provided solutions of polarized in- 
tensity from a slab of Faraday-rotating medium mixed 
with relativistic electrons that produce synchrotron emis- 
sion. The polarized intensity is expressed as a com- 
plex quantity V — Q + jU fSpangler 1982), where 
Q and U are Stokes parameters that describe the lin- 
ear polarization. Follo wing the conventions adopted 
bv ISokoloff eTall dHH), V, Q, and U are normalized 
by the total intensity of the synchrotron radiation, so 
that the modulus of V is the fractional polarization. 
The polarization angle of the complex polarization is 
* = i arctan(C//(3). 

At higher frequencies, a contribution from thermal 
emission to the total intensity reduces the observed frac- 
tional polarization by a factor (1 — fth), where fth is 
the fraction of the total intensity generated by ther- 
mal emission. We includ e thermal emission by setting 
fth = 0.23 at 4.8 GHz ([Condon fc YinI [l990t ICondod 
Il992t iNiklas et al.|[l997[ ). Determination of fth in indi- 
vidual galaxies is difficult because of the need to separate 
thermal and non-thermal emission in the radio spectrum. 
The rms scatter in fth from galaxy to gala xy is believed 
to be at most a factor of 2 ( Condonlll992D . This trans- 
lates into a spread of '--^ 25% in {I — fth) from the adopted 
value of 0.77. 

In this paper we adopt a solution for the polarized 
intensity from a unif orm slab of finite th ickness along 
the line of sight from ISokoloff et all (|1998f ). For a given 
inclination i, the component of the regular magnetic field 
along the line of sight and perpendicular to the line 
of sight B± can be evaluated. The uniformity of the 
magnetic field is expressed as the ratio of the isotropic 
random magnetic field strength to the regular magnetic 
field strength, defined as fs — ctbVS/B, where ctb is the 
rms strength of the random component of the magnetic 
field in one direction. 

The complex polarization can be written as the prod- 
uct of a wavelength-independent part that describes the 
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TABLE 4 
Sensitivity to model parameters'^ 



Inclination 


B 


/b 


Tie 


f^ 


^turb 


2h 


i = 40° 
i = 60° 
i = 75° 


0.007 
0.17 
0.50 


1.01 
1.09 
0.83 


0.006 
0.17 
0.51 


0.034 
0.12 
0.45 


0.042 
0.057 
0.055 


0.035 
0.006 
0.45 



^ Numbers listed are normalized contributions to the total scat- 
ter in a Monte Carlo simulation with 10000 realizations per in- 
clination per parameter, where all parameters were drawn in- 
dependently from a uniform distribution. Parameters with a 
small contribution to the total scatter arc relatively unimpor- 
tant. Parameter ranges adopted in this simulation: 1 < B < 10 
fiG, 0.7 < /b < 4.0, 0.01 < ne < 0.1 cm"^, 0.1 < fi < 1.0, 
10 < Iturb < 100 pc, 200 <2h < 2000 pc. 

polarization if Faraday rotation effects can be neglected, 
and a wavelength-dependent part that includes Faraday 
rotation effects. The wavelength-independent part of the 
complex polarization is 



{-Po) = 



Bi 



2a| 



exp [2j ( i TT + arctan(By /B^ ) )] , ( 1 ) 



where Pi « 0.75 is the intrinsic maximum polarization for 
synchrotron emission in a regular magnetic field, and B^ 
and By are the components of the regular magnetic field 
along two orthogonal axes in the plane of the sky. Equa- 
tion [T] assumes a constant spectral index of — 1 for the 
synchrotron emission of a galaxy. We define the x — axis 
along the apparent major axis of the disk. The brackets 
(•) indicate that the emission is averaged over a volume in 
the disk, that is sufficiently narrow to keep the regular 
magnetic field constant inside the volume, but extends 
along the line of sight through the disk. This is a rea- 
sonable approximation except for very high inclinations 
i > 80° , where the azimuthal field must change along the 
line of sight, but only if the line of sight is near the major 
axis. 

Equation [T] includes the effect of decreased polariza- 
tion because of the presence of a random component of 
the magnetic field in the plane of the sky. iSokoloff et al.l 
(|T998) referred to this as wavelength-independent depo- 
larization. The term beam depolarization is sometimes 
used in this situation, but this term is used also in case 
of wavelength-dependent Faraday effects. The plane of 
polarization is perpendicular to the component of the 
regular magnetic field projected on the plane of the sky. 

The observed polarization {V) differs from (Vq) be- 
cause the polarized synchrotron emission generated in 
a small volume element in the disk is Faraday rotated 
as it travels through the disk to the observer. The ef- 
fects of Faraday rotation can be written as a wavelength- 
dependent factor in the complex polarizat ion, t h at in - 
clu des an integral along the line of sight. iBurnI (|1966D 
and lSokoloff et al.] (|1998l ) evaluated this integral analyti- 
cally for a uniform slab. The observed polarized intensity 
from a small part of the disk is 



1 - exp(-5) 



(2) 



2jX^TZ is a complex number that 



where S — 2X'^(t^^ 
includes rotation measure fluctuations ctrm in the real 
part and the Faraday depth TZ in the imaginary part. 

The real part of S contains fluctuations in rotation 
measure along the line of sight through the factor ctrm, 




1 -1 

Relative Right Ascension 

Fig. 4. — Representation of polarized intensity in an annulus of 
a spiral galaxy disk for i = 40° (a), and i = 60° (b). Thin lines 
represent the component of the magnetic field perpendicular to the 
line of sight, B^. The pitch angle of the field is 15° in this exam- 
ple. Thick lines represent the polarization E-vectors, with length 
proportional to the polarized intensity. The line-of-sight compo- 
nent of the magnetic field By is indicated by open circles (towards 
the observer) and filled circles (away from the observer) with size 
proportional to _B|| . The integrated fractional polarization in (a) is 
EIo = 6.1%, and in (b) EIo = 13.1%, including a thermal fraction 
fffi = 0.23 at 4.8 GHz. The separation of the positions shown in 
this^ figure does not represent the step size of the integration. 

which represents the rms amplitude of rotation measure 
fluctuations after integration along the line of sight. The 
plane of polarization makes a random walk because of 
fluctuations in density and magnetic field. When polar- 
ized emission is integrated over different lines of sight, 
or different depths inside the source, depolarization of 
the emission results. This process is called internal Fara- 
day dispersion. Evaluation of ctrm requires assumptions 
about the correlation sc ale of fluctuations in density and 
magnetic field strength. ISokoloff et al.l ()1998[ ) discuss the 
difficulties estimating a correlation scale length, and as- 
sume that the scale length is the same for fluctuations 
in density, magnetic field, and rotation measure. The 
strength of the magnetic field in a turbulent cell defined 
by the correlation length is ctb- Rotation measure fluc- 
tuations also depend on the path length L through the 
disk with thickness 2h, observed at inclination i, given 
by i = 2h/ cosi. The line-of-sight flUing factor of the 
Faraday- rotating medium is fi. Rot ation measure fluctu - 
ations arc then calculated foUowing lSokoloff et al.l (|1998l ) 



Ctrm — O-SfrieCTB^turb 



\lu,rh cosi/ 



(3) 



The imaginary part of S describes rotation of the 
plane of polarization by the line-of-sight component of 
the regular magnetic fleld. The Faraday depth TZ — 
nefiB\]2h/ cosi is the product of the mean electron den- 
sity (riefi), where is the electron density in a turbulent 
cell, the line-of-sight component of the regular mag- 
netic flcld, and h the path length through the disk. Rota- 
tion of the plane of polarization also depolarizes emission 
integrated along the line of sight as emission from differ- 
ent depths of the source is rotated by different amounts. 
This mechanism for depolarization is called differential 
Faraday rotation. It is sometimes called depth depolar- 
ization, but it should be understood that there is a differ- 
ent depolarizing effect in the form of Faraday dispersion 
that is also the result of integrating emission along the 
line of sight. 

3.2. Integrated polarized emission 
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TABLE 5 
Model parameters for Figure [T] 



Model 


B 


/b 




/» 


^turb 


2h 


fth 


n-niax 


nmcd 








(cm^'') 




(pc) 


(pc) 




(%) 


(%) 


1" 


5.0 


1.0 


0.03 


0.5 


50 


1000 


0.23 


16.4 


10.2 




10.0 


1.0 


0.03 


0.5 


50 


1000 


0.23 


13.7 


8.2 


T 


5.0 


1.0 


0.09 


0.5 


50 


1000 


0.23 


12.0 


7.0 



4<i 


5.0 


2.0 


0.03 


0.5 


50 


1000 


0.23 


6.8 


4.7 


5" 


10.0 


2.0 


0.03 


0.5 


50 


1000 


0.23 


5.6 


3.7 


6" 


5.0 


2.0 


0.09 


0.5 


50 


1000 


0.23 


4.8 


3.1 


ja 


5.0 


3.0 


0.03 


0.5 


50 


1000 


0.23 


3.4 


2.4 


8" 


10.0 


3.0 


0.03 


0.5 


50 


1000 


0.23 


2.7 


1.8 


9'= 


5.0 


3.0 


0.09 


0.5 


50 


1000 


0.23 


2.2 


1.4 



^ Models indicated by a solid curve in Figures [1] and [If Models indicated by a 
dashed curve in Figui'es^and[5f Models indicated by a dotted curve in FigurcsfTI 
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Equation [5] provides the complex polarized intensity 
(V) for every line of sight, from which we obtain the 
local Stokes Q and U intensities, normalized by the local 
total intensity that is presumed uniform in the current 
models. The model Stokes Q and U intensities can be 
integrated over the disk to predict the model polarized 
flux density, fractional polarization IIo, and the angle 0b 
of an unresolved spiral galaxy. The integration over the 
disk reduces to a one-dimensional numerical integration 
over azimuthal angle for the symmetric models under 
consideration. 

Assuming values for the various model parameters, the 
integration over the disk can be performed as a function 
of inclination. The integrated polarization is more sen- 
sitive to certain model parameters, and the influence of 
model parameters can change with inclination. The rel- 
ative importance of the parameters in our model was 
investigated through a Monte Carlo experiment. Each 
parameter was varied separately, and the resulting range 
of the integrated polarization was divided by the range 
obtained when all parameters were varied independently 
at the same time. The spread in integrated fractional po- 
larization was taken as the difference between the upper 
quartile and the lower quartile. As the underlying distri- 
butions are strongly skewed, the quartiles provide a bet- 
ter indication of the spread in the distribution than the 
standard deviation. The normalization with the "total" 
scatter allows us to compare the importance of parame- 
ters between the three inclinations. It should be noted 
that the non-linear character of the model, and the op- 
posing effects of certain parameters can make the ratio 
of the normalized quartile ranges larger than unity. 

Table [5] shows for three values of the inclination the re- 
sult of 10000 realizations of each parameter, drawn from 
a uniform distribution. The magnetic field strength, elec- 
tron density, filling factor, turbulent length scale, and the 
vertical thickness of the disk were varied by an order of 
magnitude, as indicated in the footnote to Table [H The 
ratio of random to regular magnetic field strength (/b) 
was varied over a smaller range (0.7 to 4.0). At low in- 
clinations, the normalized quartile range for variation of 
/b is 1.01, while the normalized quartile ranges of all 
other parameters are less than 0.1. Clearly, the ratio 
of random to regular magnetic field dominates the frac- 
tional polarization of the model for small inclinations. 



At i = 60°, /b is still the most important parameter, 
but other parameters such as the regular magnetic field 
strength and the electron density become more impor- 
tant. At i — 75°, /b is still important, but the parame- 
ters that define TZ = nefiB\]2h/ cosi are almost equally 
important. The relative importance of rig, /i, B, and 
2h is similar at i = 75°, because they contribute equally 
to 7?., although subtle differences occur because of the 
dependence on ctrm- The turbulent scale length Iturb is 
the parameter in our model that is least constrained by 
observations. Table 2] shows that its effect on our results 
is insubstantial. 

The simplest model is a disk with uniform parameters 
throughout that can be represented by a set of concentric 
annuli, all with the same integrated polarization proper- 
ties. The model makes no assumptions about the mag- 
netic field structure in the center of the disk. This is not 
a significant simplification for most spiral galaxies, where 
most of the polarized emission originates in the disk. Fig- 
ure |4] shows a graphic representation of polarized inten- 
sity at a representative subset of locations (a ring) in the 
disk for model 1 in Table [H] at inclinations i = 40° and 
i = 60°. Thin lines indicate the magnitude and direction 
of the magnetic field component in the plane of the sky, 
B±. Thick lines represent polarization vectors. Open 
and filled circles represent the component of the regular 
magnetic field along the line of sight, B\^ . 

At low inclination (Figure[4^), the path length through 
the disk is relatively small and the line-of-sight compo- 
nent of the regular field is also small. The polarized in- 
tensity in Figure |4^ is fairly uniform across the disk and 
relatively high. At higher inclination, the path length 
through the disk increases, and the line-of-sight compo- 
nent of the regular magnetic field increases, in particular 
near the major axis. The polarized intensity is smaller 
everywhere in the disk, but most strongly near the major 
axis. Here, the line-of-sight component of the azimuthal 
regular field is largest, so Faraday rotation is strongest 
and the component B± (Equation[T]) is smallest. Faraday 
depolarization is therefore strongest near the major axis, 
and the intrinsic polarization of the emission is smaller 
than elsewhere in the disk. Faraday rotation of the plane 
of polarization is also stronger near the major axis in Fig- 
ure g] 

The model Stokes Q and U intensities, calculated by 
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solving Equation [2] for every line of sight through the 
disk, are integrated over the disk to predict the polarized 
intensity of an unresolved spiral galaxy. The integration 
over the disk reduces to a one-dimensional numerical in- 
tegration over azimuthal angle for the symmetric models 
under consideration. The model in Figure has in- 
tegrated fractional polarization XIq = 7.9%, while the 
model in Figure [4}3 has integrated fractional polarization 
Ho ~ 17.0%. Despite the stronger depolarization any- 
where in the disk at inclination i = 60°, the integrated 
polarization is stronger at z = 60° than at i = 40°, be- 
cause of the smaller degree of symmetry of the polarized 
emission. 

The model predictions for IIo and are independent 
of the spiral pitch angle of the magnetic field. For a disk 
with constant thickness and uniform properties, the path 
length through the disk depends only on inclination, and 
the polarized intensity depends only on the angle of the 
regular component of the magnetic field with the line of 
sight. The effect of a pitch angle of the magnetic field, is 
to change the azimuthal angle in the disk where the line 
of sight makes a certain angle with the regular magnetic 
field. The integral over all azimuthal angles in the disk is 
not affected by a translation in azimuthal angle resulting 
from the pitch angle of the magnetic field. 

Table [5] gives values for the parameters in the models 
that are compared with th e data . Mode l 7 corresponds 
with the values adopted bv lBeckl ()2007b[ ) for NGC 6946. 
The filling factor of ionized gas ji and the thickness of the 
disk 2h were kept constant, because their effect on the 
models is similar to that of the electron density and the 
magnetic field. Model curves of fractional polarization 
as a function of inclination are shown in Figure [TJ Three 
families of models are shown (Table [5]), with /b = 1, 2, 
and 3. For each value of /b, model curves are shown for 
Tie ~ 0.03 cm"'', B = h /iG (solid curves), Ue — 0.03 
cm^'^, B = 10/iG (dashed curves), and Ue — 0.09 cm~'', 
B = 5 /iG (dotted curves). 

At low inclination, the models are mainly distinguished 
by the value of /b. The line-of-sight component of the 
regular magnetic field is small at low inclination, and 
the main depolarization mechanisms are wavelength in- 
dependent depolarization and Faraday dispersion. At in- 
clination i > 50° , differential Faraday rotation along the 
line of sight becomes increasingly important, and signif- 
icant effects of the strength of the regular magnetic field 
and the density of thermal electrons are found. 

This effect already occurs at low inclination, where the 
line-of-sight component of the magnetic field near the 
major axis increases proportionally to i, but the path 
length through the disk ~ 1/ cos i is constant to first or- 
der in i. The result is that models that include Faraday 
depolarization predict a higher degree of polarization at 
low inclinations th an models that do not include Fara- 
day depolarization (jStil et al.ll2007f ). The present models 
also predict a higher fractional polarization at low incli- 
nation if the electron density or the strength of the 
regular magnetic field is increased, if the ratio of random 
magnetic field to regular magnetic field strength is not 
too high. 

The models more or less cover the same area of Fig- 
ure [T^ as the data. This result supports the idea that 
the observed integrated polarization of spiral galaxies is 
related to the large-scale field in the disk. The most no- 
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Ho (%) 

Fig. 5. — Probability density functions (PDF) of Ho for 
randomly-oriented spiral galaxies derived from the models in Fig- 
ure [T] The continuous, dashed, and dotted curves correspond with 
models 1 - 9 in Table [S] in order of decreasing maximum Ilmax ■ 
The curves are normalized such that the integral of the PDF is 1. 
Median values Ilnied of these distributions are listed in Table \5\ 

table difference between the data and the models is the 
observed low Hq < 2% of some galaxies with intermediate 
inclination, 40° < i < 60°. It is difficult to change any 
single parameter in the model to obtain a fractional po- 
larization n < 2% in this inclination range, unless higher 
values of /b are assumed. 

3.3. Probability distribution of Hq 

The results in the previous sections suggest that the 
integrated 4.8 GHz emission of spiral galaxies can show 
substantial polarization. Deep polarization surveys can 
detect distant unresolved spiral galaxies. How many 
spiral galaxies can be detected, and the interpretation 
of unresolved galaxies detected in polarization, is de- 
termined by the probability density function (PDF) of 
the fractional polarization of spiral galaxies. The wave- 
length dependence of integrated polarization is also best 
summarized by means of the PDF of fractional polariza- 
tion. In this section we discuss the expected shape of 
the PDF of fractional polarization for unresolved spiral 
galaxies at 4.8 GHz and 1.4 GHz. We also discuss the 
PDF of Ho in the special case of low-inclination galax- 
ies {i < 30°), where the integrated polarization depends 
mostly on Faraday dispersion, but the actual value of 
the inclination cannot be determined reliably from the 
optical axial ratio. 

The models for the integrated fractional polarization 
as a function of inclination can be used to derive a prob- 
ability density function (PDF) for Hq of unresolved, ran- 
domly oriented spiral galaxies. This PDF can be used to 
predict the number of spiral galaxies detectable in deep 
polarization surveys based o n population rn odels of the 
total intensity source counts (jStil et al.ll2007f ). The shape 
of the PDF can be modeled as a function of rest-frame 



10 



Stil et al. 



2000 - 



CO 

*x 
o 

CO 



1000 



X5 







J I i_i I I i_i I I L.I I I I L 



(a) 



T|I|TTnnnTTWT 



2000 - 



CO 

x 
o 

o 

CD 



^ 1000 



J3 

13 







_l I I I I I I I l_ 



(b) 



Tl]lTn7niT|rTT>^^ , ,, 







10 15 20 25 







10 



Ho (%) 



Ho (%) 



wavele ngth and fitted to observations. iTavlor et al.l 
(|2007| ) used Monte-Carlo simulations of noise statistics 
and observational selection effects, and derived the shape 
of the PDF for XIq for faint polarized sources from a max- 
imum likelihood fit to the data. A similar approach can 
be used in the future to compare the PDF of Ho for 
galaxies at high rcdshift with the PDF for local galaxies. 

The PDF for XIq follows from the probability P{i) to 
observe a galaxy at inclination between i and i + di, 



P{i)di = siii{i)di (0° < i < 90°). 



(4) 



Fig. 6. — (a) Simulated histogram of fractional polarization of 20,000 randomly oriented model spiral galaxies at 4.8 GHz. Ranges of 
model parameters used for this simulation arc given in Table |4] The fractional polarization for the same galaxies at 1.4 GHz is shown in 
(b). The vertical axis in (b) is truncated (lowest bin only). 

Qualitatively, the shapes of the PDFs in Figure [5] do 
not depend strongly on the details of the models. The 
PDFs are derived from the curves in Figure [TJ The 
shapes of these curves can be understood qualitatively 
in terms of two competing effects. If the inclination is 
not too large, the integrated polarization increases with 
inclination because of a larger net magnetic field in the 
plane of the sky. At higher inclinations, Faraday de- 
polarization dominates and the integrated polarization 
decreases. A broad maximum in XIq must occur at in- 
termediate inclination, that gives rise to the peak of the 
PDF near the maximum of XIq. 

Another feature of the PDFs in Figure [5] is that they 
converge to a constant value for small Hq. The distri- 
bution is dominated here by galaxies at low inclination, 
with a smaller contribution of highly inclined galaxies 
with strong Faraday depolarization. Small fluctuations 
in the curves at small Ho are the result of occurrences 
where the Faraday rotation of the disk at the major axis 
is TZX^ = mr {n = 1, 2,3,.. .) at high inclination. These 
fluctuations are small, and would not be noticed in ob- 
servations. 

The plateau in the distributions in Figure [5] can be 
understood based on a simple analytic argument. For 
inclination i < 30°, the model curves in Figure [T] can 
be approximated with the relation Hq = ai'^, with the 
inclination i in radians and a a constant that depends on 
the model. The value of a may depend on complicated 
physics, but the only assumption here is that the increase 
of Hq for small i can approximated by a parabola, for 



The probability to find a fractional polarization between 
Hq and XIq -I- AXXq is equal to the probability to find a 
galaxy with inclination in one of up to two inclination 



ranges 



to Ir, 



that yield the same fractional polar- 



ization. Therefore, 



N 



F(no)Ano =^(cos 



^min.fc COS Imax.k ) 



(5) 



k=l 



In most cases, N = 2, except at the maximum of XIq, 
where A'" = 1. 

The PDFs of Ho for models 1 - 9 in Table [5] are 
shown in FigurelH The median inclination for randomly- 
oriented galaxies is «mcd = 60°. Nearly half (47%) of the 
galaxies in a randomly oriented set have inclination be- 
tween 50° and 80°, for which XXq is near the maximum 
of the curves in Figure [T^. Consequently, all models pre- 
dict a PDF that is sharply peaked toward the maximum 
possible value of XXq. 
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fixed parameters of the model. The probability to find a 
fractional polarization between Ho and Ho + cHIq is equal 
to the probability to find a galaxy with inclination in 
the interval ii ^ i ^ 12 < 30°, with Eq = ai^i^,!, and 
(Ho + rfllo) = a«max,i- Equation [5] reduces to 

1 / sin. z \ 

P(no)dno = COSZniina-COSlmax,! = 2a V^~/ 

The last step used that (imax,i~*min,i) is the infinitesimal 
increment di = dllf)/{2ai), and that ^(imin,i + *max.i) = 
i. For i < 30° the factor sini/i « 1 to an accuracy 
better than 5%, so ^'(no) is a constant that depends only 
on a. The models presented in Section [3] suggest that 
a depends mainly on the ratio of the random magnetic 
field to the regular magnetic field. This is why models 
1, 2, and 3 converge to the same value for small XIq, 
and models 4, 5, and 6 converge to a larger value. The 
differences between models 7, 8, and 9 arise from the 
more significant contributions from galaxies with a high 
inclination. These models also converge to a common 
value if only galaxies with low inclination are considered. 

The distribution of Ho of an observed sample of spiral 
galaxies would be a superposition of a continuous set of 
model distributions similar to the nine curves in Figure^ 
A simulated PDF for 20,000 galaxies with model param- 
eters varied over the ranges given in Table 3] is shown 
in Figure [6] At 4.8 GHz (Figure [6^), the distribution 
of Ho is quite broad, with a peak between 1% and 2% 
polarization. The location of the peak and the depth of 
the minimum near Hq = 0% depends on the number of 
galaxies with strong Faraday depolarization. At 1.4 GHz 
(Figureinb), Faraday depolarization is much stronger and 
the distribution is sharply peaked towards Hq = 0%. For 
the distributions in Figure [6] we find that 78% of galax- 
ies has Ho > 1% at 4.8 GHz, while at 1.4 GHz 17% has 
Ho > 1%. The broad tail of the distribution in Figure^Js 
contains galaxies with low to moderate inclination. The 
median inclination for galaxies with Hq > 5% at 1.4 GHz 
is 41°. 

Integrated polarimetry of an unbiased sample of a few 
hundred spiral galaxies would allow a statistical analy- 
sis for the Ho distributions of spiral galaxies, including 
subsets by inclination, morphological type and star for- 
mation rate. The Ho distribution can be retrieved from 
the data in the presenc e of noise through m aximum like- 
lihood fits presented bv lTavlor et al.l ()2007l ). The skewed 
shape on the PDF in Figure [S^, suggests that the best 
theoretical function to represent the PDF of spiral galax- 
ies would be the lowest order onti-symmetric deviation 
from a Gaussian that shifts the peak of the distribution to 
a finite value of Hp. I n term s of the Gauss-Hermite series 
used in lTavlor et al.l ()2007l), the lowest order term w ould 
include coefficient /13 (|Van der Marel fc Franxl[l993l ). 

4. DISCUSSION 

4.1. Diagnostics from integrated polarization 

We presented observations of the integrated polariza- 
tion of nearby spiral galaxies to study the polarization 
properties of spiral galaxies as unresolved radio sources. 
It is shown that unresolved spiral galaxies are polarized 
radio sources with fractional polarization up to ^ 20% 
at 4.8 GHz. This high degree of fractional polarization 
arises from the regular azimuthal component of the mag- 



netic field in the spiral galaxies observed, that is pro- 
jected into a component in the plane of the sky that is 
predominantly oriented along the apparent major axis of 
the disk seen at inclination i. The projected magnetic 
field direction as implied by the plane of polarization, 
rotated by 90° , is well correlated with the position angle 
of the major axis for galaxies observed at higher incli- 
nation. The highest integrated polarization is expected 
for galaxies at intermediate inclination. For a sample of 
randomly oriented spiral galaxies, 47% has an inclination 
between 50° and 80°. Galaxies in this inclination range 
are under-represented in the present sample because the 
data used in this paper were obtained for polarization 
studies that preferred low-inclination or high-inclination 
galaxies. The number of galaxies with a high fractional 
polarization in Figure [T] may therefore be under repre- 
sented compared with an unbiased sample. 

Three of the eight Virgo spirals in our data set have 
integrated fractional polarization Hq ^ 8%. The range 
of Ho of the Virgo spirals is similar to that of the nearby 
"field" spirals. We do not find evidence in the present 
data that the integrated polarization of spiral galaxies 
in a cluster environment would be different from field 
galaxies. Although cluster galaxies are distorted by inter- 
actions with other galaxies and the intraclustcr m edium 
(.Wez gowiec et al. 2007; Chyzy_et al. 2007; Vollmcr eTall 
|2007( ). the Virgo galaxies in our sample cannot be distin- 
guished from field galaxies in their integrated polariza- 
tion. This result should be verified with a larger, com- 
plete sample of cluster galaxies. 

The barred galaxies in Table [3] show that a high in- 
tegrated polarization also occurs for barred galaxies. In 
some of these barred galaxies, polarized emission from 
the bar is negligible compared to polarized emission from 
the disk. In other barred spirals, polarized emission from 
the bar region dominates, and the integrated polarization 
depends on the orientation of the bar, as well as the incli- 
nation of the galaxy. The effect of a bar is not included 
in the present models, so the existence of significant de- 
partures from the models is not unexpected. Non-axially 
symmetric models of the integrated polarization of spiral 
galaxies are deferred to a later paper. 

Figure [3] suggests a difference between luminous and 
less luminous spiral galaxies. The three sub samples sep- 
arately also show a smaller integrated fractional polar- 
ization for luminous spiral galaxies. Face-on galaxies and 
depolarized edge-on galaxies should exist at every lumi- 
nosity, although they form a minority in a set of ran- 
domly oriented spiral galaxies. The presence of some 
galaxies with low Hq at every luminosity in Figure [3] is 
therefore not surprising. The luminous galaxies in the 
sample cover the same range in inclination as the galax- 
ies with lower luminosity. This excludes selection effects 
in inclination as a possible origin of the correlation. No 
selection effect could be identified that would exclude lu- 
minous galaxies with a high Hq from the sample. 

The lack of luminous galaxies with a high Hq gains sig- 
nificance in view of the results from Section 13.31 Given 
certain physical conditions, the inclination dependence 
of Ho favors high values of Ho. The lower polarization 
of luminous spiral galaxies therefore indicates different 
conditions in these galaxies. A higher rate of star for- 
mation increases the random component of the magnetic 
field through the effect of supernova explosions and stcl- 
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lar winds. A stronger total magnetic field and lower po- 
larization are found locally in regions of galaxies wher e 
the star formation rate is high (|Beckll2005l : IChvzvll2008D . 
The correlation in Figure [Sja may be explained similarly 
in terms of a generally higher star formation rate. The 
correlation of polarization with luminosity is an impor- 
tant factor for predictions of the number of spiral galax- 
ies that will be detected in deep polarization surveys. 
Larger samples of spiral galaxies are required to confirm 
this result, and to allow more detailed modeling. 

Simple axially symmetric models for the integrated po- 
larization of spiral galaxies successfully describe the dis- 
tribution of observed galaxies in the Hq - i diagram (Fig- 
ure [TJi). The results in Table [1 show that for i < 50° 
the model is most sensitive to the parameter /b- This 
is so because wavelength-dependent depolarization such 
as Faraday dispersion and differential Faraday rotation 
are small compared with wavelength-independent depo- 
larization in the models at 4.8 GHz with a small incli- 
nation. At higher inclination and at lower frequencies, 
wavelength-dependent depolarization is more important. 
A complete sample of spiral galaxies observed at 4.8 GHz 
may be divided into sub samples by inclination to sep- 
arate the effects of different physical parameters. We 
discuss three inclination ranges. 

Galaxies with inclination i < 30° (13.4% of an un- 
biased sample) would display low fractional polarization, 
provided the galaxy is axially symmetric. Large-scale de- 
partures from axial symmetry likely result in a higher de- 
gree of polarization for these low-inclination galaxies. A 
statistical study of the integrated polarization of face-on 
spiral galaxies will therefore provide constraints on the 
large-scale asymmetry of polarized emission and there- 
fore magnetic fields in spiral galaxies. For a sample of 
face-on barred galaxies, the orientation of the plane of 
polarization may be correlated with the position angle of 
the bar derived from optical images. 

At inclinations 30° < i < 50° (22.3% of an unbiased 
sample), the axially symmetric models are most sensi- 
tive to the ratio of the random magnetic field to the reg- 
ular magnetic field. This inclination range seems most 
promising to provide constraints of this parameter for 
a large sample of galaxies. Integrated polarimetry at 
more than one frequency may be combined to separate 
wavelength-dependent depolarization from wavelength- 
independent depolarization. More detailed models of the 
effect of non-axially symmetric structure such as spiral 
arms are required. 

At inclinations i > 50° (64.3% of an unbiased sam- 
ple), the longer path length and larger line-of-sight com- 
ponents of the regular magnetic field greatly enhance 
wavelength-dependent Faraday depolarization and Fara- 
day rotation. This is confirmed by the results in Tabled 
where the contribution of B, rig, fi, and 2h to the total 
scatter all rise to approximately equal amounts at incli- 
nation 75°. The ratio of random to regular magnetic field 
is still an important factor, but this can be constrained 
by results obtained from galaxies at lower inclination. 
Multi-frequency observations would be valuable because 
the models suggest that the strongest frequency depen- 
dence of depolarization of the integrated emission occurs 
for inclinations more than ~ 60°. 

4.2. Spiral galaxies as polarized background sources 



Are spiral galaxies suitable background sources for ro- 
tation measure studies? Internal Faraday rotation in dis- 
tant spiral galaxies likely amounts to rotation measures 
of a few hundred rad m~^. However, the correlation be- 
tween 0B and 0opt in Figure[2^, and the models presented 
in Section [3] suggest that the plane of polarization is re- 
lated to the orientation of the major axis of the opti- 
cal disk. The polarization angle of an unresolved spiral 
galaxy is determined by the orientation of its major axis 
through the projection of the azimuthal magnetic field in 
the galactic disk on the plane of the sky. This is in part a 
result of the symmetry of the galaxies, as polarized emis- 
sion from areas of the disk opposite to each other with 
respect to the center of the galaxy is Faraday rotated by 
equal but opposite angles. Depolarization resulting from 
differential Faraday rotation also mitigates the effect of 
Faraday rotation on the position angle of the integrated 
emission. The correlation in Figure [2^ confirms that the 
plane of polarization is related to the orientation of the 
major axis in real galaxies. 

The important implication is that the plane of polar- 
ization of a spiral galaxy that is unresolved by a radio 
telescope, is independent of wavelength. The orienta- 
tion of the plane of polarization can be determined from 
the position angle of the optical major axis, although 
with a 180° ambiguity. As the plane of polarization does 
not depend on wavelength, we find the surprising result 
that unresolved spiral galaxies behave as idealized polar- 
ized background sources without apparent internal Fara- 
day rotation. This result follows from axial symmetry. 
The 8.4 GHz polarization data in Table [2] seem to con- 
firm that the polarization angle of the integrated emis- 
sion is independent of wavelength. More observations 
are required to confirm this mainly theoretical result for 
more galaxies and at lower frequencies. 

The intrinsic position angle of the plane of polarization 
of an unresolved spiral galaxy can be estimated from an 
optical image. It does not depend on internal Faraday 
rotation or on the spiral pitch angle of the magnetic field 
(Section [3|) . In principle, this property allows the de- 
tection of possible Faraday rotation by an intergalactic 
medium, using a large sample of spiral galaxies. The cor- 
relation between 9b and 6'opt (Figure [2]) can be made for 
a large number of galaxies, subdivided by distance and 
location in the sky. If large-scale magnetic fields exist in 
the nearby intergalactic medium, the correlation between 
and ^opt would show an offset for more distant galaxies 
in a certain direction. The rotation measure sensitivity 
of such an experiment would depend primarily on the 
sample size of spiral galaxies with integrated polarime- 
try, and on the longest wavelength at which unresolved 
spiral galaxies remain polarized sources. 

5. SUMMARY AND CONCLUSIONS 

This paper discusses the integrated polarization prop- 
erties of nearby galaxies from archival data at 4.8 GHz, 
and models for the integrated polarized emission of a 
spiral galaxy. We find that 

1. Spiral galaxies can show substantially polarized (up 
to ~20% at 4.8 GHz) integrated radio emission. The 
polarization angle of the integrated emission of spirals 
without a (strong) bar is aligned with the apparent mi- 
nor axis of the disk. Barred spirals can also be polarized 
sources, but do not show a clear correlation of polariza- 
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tion angle with the orientation of the major axis. 

2. The highest degree of polarization, and the largest 
spread in fractional polarization, are found for spiral 
galaxies with intermediate to high inclination (i > 50°). 

3. A sub-sample of nine spiral galaxies in the Virgo 
cluster does not show a significant difference in inte- 
grated polarization properties from the field galaxies in 
our sample. This suggests that physical distortions of 
these galaxies are not so strong that they have a notice- 
able effect on the integrated polarization. Confirmation 
by a larger, complete sample is required. 

4. The fractional polarization of integrated emission is 
smaller for galaxies with a higher radio luminosity. This 
result suggests a difference in properties of the magnetic 
field and/or Faraday depolarization in luminous spiral 
galaxies. 

5. Our models for the integrated polarization of axially 
symmetric spiral galaxies cover the same area in the IIo- 
i diagram (Figure [T]) as the observed galaxies without a 
strong bar. The models indicate that the main cause of 
variation in the fractional polarization of galaxies with 
inclination less than 50° is the ratio of the random to 
the regular magnetic field component in the disk. 

6. Our models provide a probability density function 
of fractional polarization of spiral galaxies as a function 
of wavelength that can be used to predict the number 
of spiral galaxies in deep polarization surveys, and to 
study the statistical properties of spiral galaxies in such 
surveys. 
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